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ABSTRACT 

In this paper we present a topological magnetic field investigation of seven two- 
ribbon flares in sigmoidal active regions observed with Hinode, STEREO, and SDO. 

We first derive the 3D coronal magnetic field structure of all regions using marginally 
unstable 3D coronal magnetic field models created with the flux rope insertion method. 

The unstable models have been shown to be a good model of the flaring magnetic 
field configurations. Regions are selected based on their pre-flare configurations along 
with the appearance and observational coverage of flare ribbons, and the model is con¬ 
strained using pre-flare features observed in extreme ultraviolet and X-ray passbands. 

We perform a topology analysis of the models by computing the squashing factor, Q 
in order to determine the locations of prominent quasi-separatrix layers (QSLs). QSLs 
from these maps are compared to flare ribbons at their full extents. We show that in 
all cases the straight segments of the two J-shaped ribbons are matched very well by 
the flux-rope-related QSLs, and the matches to the hooked segments are less consistent 
but still good for most cases. In addition, we show that these QSLs overlay ridges in 
the electric current density maps. This study is the largest sample of regions with QSLs 
derived from 3D coronal magnetic field models, and it shows that the magnetofrictional 
modeling technique that we employ gives a very good representation of flaring regions, 
with the power to predict flare ribbon locations in the event of a flare following the time 
of the model. 


Subject headings: Summagnetic fields — Sun: sigmoid — Sun: flares — Sun: X-rays 
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Introduction 


Solar flares, often accompanied by coronal mass ejections (CMEs), are the most powerful 
phenomena originating from the variable sun. The energy powering these events is thought to 
come from magnetic reconnection, a process that converts fre e magnetic energy into radiation, 


energetic particle acceleration, and kinetic energy of plasma ([Forbes et ah 2008 ). Various ob¬ 


served flare and CME-related features can be understood in the context of the standard 2D flare 


mo del, the basics which wer e laid out by ICarmichaell (j 1964 1 , ISturrockl (|1966l l , iHiravamal (119741. 


andlKonn fc Pneumanl (119781. an d further extended based on Yohkoh observations bv IShibata et aL 
(Jl995l l and It sunetnl (|l996l . Il997l l. In this model, flare ribbons correspond to the chromospheric 


footpoints of newly-reconnected field lines, and flare loops represent reconnected field lines cool¬ 
ing beneath the reconnection site at the base of the escaping flux rope. While the 2D model has 
successfully explained multiple related phenomena, flares are intrinsically 3D events. The physical 
processes underlying flares, particularl y magnetic reconnection, may have pro perties in 3D that 
significantly deviate from the 2D view (jAulanier et al .1120061 : iJanvier et al.ll2013l l. For example, the 
shape of the ribbons and the existence of transient coronal holes (dimmings), associated with the 
footprints of the escaping flux rope, are only explained by the 3D models. 

It is therefore necessary to better understand the 3D properties of ultraviolet (UV) flare rib¬ 
bons because of their complex morphology and dynamics offer a particularly straightforward way 
to study the three-dimensionality of flares. UV flare ribbons are produced by heat conduction and 
energetic pa rticles streaming dow n the flare loops, which interact with the denser layers of the solar 
atmosphere ([Fletcher et al,J[201l])- Flare ribbons arise in the process of chromospheric evaporation, 
whereby chromospheric heating from the corona is suddenly increased and overcomes the cooling 
rate of the chro mosp h eric material , dramatically increasing the density and temperature of upflow- 
ing plasma (lAntonucci et al.lll984T). This pr o cess c an be aided by particle precipitation from the 
reconnection site, as discussed in llteid et al.1 ({2912 ). Because heat conduction fronts and particles 
flow down magnetic field lines, the shape, separ ation, and size of fl are ribbons are intimately related 
to the magnetic structure under the flare site (lY urchysh yn et al. 2008 ). 


Although chromospheric evaporation is included in the 2D standard model, both for producing 
ribbons and populating the flare loops with dense plasma, a 2D view cannot facilitate predictions 
about the shape of flare ribbons. When extending the standard model concepts to 3D, the shape 
of the ribbons are supposed to be determined by the footpoint locations of reconnected field lines, 
along which energy is flowing down. The shape and dynamics of the flare ribbons are thus tightly 
linked with the 3D magnetic configuration and dynamics during a flare. Ribbons can therefore 
be used to identify properties of the reco nnection site, to gai n insi ght into the particular flare 
configuration, and to test 3D flare models ([Gorbachev fe Somov 19881 ). 


Magnetic topology aims to studying and identify specific locations in the 3D magnetic field 
where the strong electric current sheets necessary to induce reconnection are likely to form. This en¬ 
ables us to determine the field lines that will be involved in reconnection and to predict along which 
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field lines energy will flow to form fla re r ibbo ns. A separatrix surface, through which the magnetic 


field connectivity is discontinuous fLongcope 2005), is such an example of a magnetic configuration 


preferential for magnetic reconnection. T he f i rst to match ri bbon s to separatrices from a potential 
field model of an active region were Gorbachev &, Somov 19881 ). Overlap between flare ribbons 


ces on the photosphere from potential and linear force-free field (LFFF) models ( 

Mandrini et al. 

1991, 

1995, 

2014: 

Demoulin et al. 

1993; 

Demoulin et al. 

1994; 

van Driel-Gesztelvi et al. 1994). A 


null point is another example of a topological feature exhibiting a distinctive shape: observations 



the Js adjacent to each other on either side of the polarity inversion line (PIL) and the hooks facing 
each other on the opposite ends. Understanding the magnetic topology associated with these two 
J-shaped ribbons can enlighten us on the processes leading to and developing during eruptions. A 
subcategory of the two J-shaped ribbon flares are those for which the active region of origin presents 
a forward or inverse S-shape in its extreme ultraviolet (EUV) or X-ray em ission for some time prior 
to the flare. These regions are commonly known as sigmoids ((Rust fe K umar 199a ). Sigmoid-region 
two J-shaped ribbon flares have been interpreted to occur in bi-polar configurations containing a 
flux rope. Evidence for the presence of flux ropes comes from magnetic modeling as well as the very 
shape of the ribbons. Non-linear force-free fiel d (NLFFF) models and extrapolations have been 
used t o reveal the presence o f a twisted flux rop e ( Schriiver et al . 200 ( 3 : Isavcheva fc van Ballegooiienl 


2009 


Savcheva et al. 2012a 3; Liu et al. 2014 ). Additionaly, the J-shape of the ribbons has been 


shown to correspond to the specific topological structures associated with a twisted flux rope: 
quasi-separatrix layers (QSLs) and hyperbolic flux tubes (HFT). 

QSLs are purely 3D topological structures that do not exist in 2D. They correspond to sub- 


experiences drastic changes ( 

Priest & Demoulin 

1995 

Demoulin et al. 

1996; 

Titov et al. 2002; 

Titov 

2007: 

Pariat & Demoulin 

20121). QSLs, as separatrices, are preferential sites for the development of 

magnetic reconnection (Demoulin et al. 1996: 

Aulanier et al. 2005, 

2006; 

Janvier et al. 20131). 

HFTs 


are lo cated at the intersection of QSLs and are a generalization o f magnetic separators (ITitov et al 


20021) th at are typical of the regions below twisted flux ropes ifritov et al 


2002 


2012al lbl: IPariat <fe Demoulin 2012 : Janvier et al. 2013 : Guo et al. 2013 : Zhao et al 


Savcheva et al 


20141). Interest¬ 


ingly, the photosp heric footprints of the QSLs surrounding the twisted flux rope pre sent the shape 
of two inverted Js ( Titov et al.ll2002l : ISavcheva et al.ll2012bl : IPariat <fe Demouliiill2012n . Magnetic re¬ 
connection occurring at these QSLs would thus naturally lead to the generation of the two J-shaped 
ribbons observed during flares. Observation s of two J-shaped flare ribbons are therefore strong evi- 


dence for the presence of a twisted flux rope ([Titov et af 


2002 


Schriiver et al. 


2011 


Savcheva et al. 
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2012bl : IZhao et al.ll2014l : iDudfk et al.112014 ). J-shaped ribbons have also been associated with ob¬ 


servations of J-shaped electric currents during flares, which is theoretically predicted b y the 3 D 
standard model since the QSLs are preferential sites for current accum ulation ( Janvier et a l. 20141 ). 
Finally, reconnection developing at 3D QSLs, slipping reconnection ((Aulanier et al. 2006|), differs 
from 2D cut-and-paste reconnection. Motions of UV brig hten in g observed during the formation of 
the J-shaped ribbons have exhibited slipping behaviors (|Dudfk et al. 2014), further confirming the 
presence of J-shaped QSLs in sigmoid flare events. 


However, most previous studies have only examined case-by-case connections between J-shaped 
ribbons and magnetic topologies, and no statistical studies have been done. A selection bias 
potentially exists for regions chosen based on their favorable topological properties, i.e. only regions 
that do present a close correlation between J-shaped ribbons and J-shaped QSLs have been studied 
and published on. In this study, we find particular topologies for a sample of well-observed events. 
We present a correspondence between the morphology of the photospheric QSLs computed from 
observationally-constrained magnetic field models of the coronal field and the flare ribbons for 
7 sigmoid region eruptions. These regions have been selected based on their X-ray and EUV 
morphology prior to the flare and on the shape of the ribbons. No topology criteria have been 
used, and we therefore reach a more statistically significant conclusion on whether or not J-shaped 
flare ribbons in sigmoid regions tend to match the distribution of QSLs. This study allows us to 
control whether a portion of J-shaped ribbons is better matched to the coronal magnetic field model 
than others. Such comparisons with observed ribbons enables us to study and control the quality 
of our magnetic field model in a new and innovative way. 


This study is the observatio nally-constrained counterpart of the simulation results of Aulanie r et al. 


next step after (jSavcheva et al 


(120121) and Ja nvier et al.1 (120131). w ho use QSLs to extend the standard flare model to 3D and the 


2012bl . ; hereafter S12). The ultimate test for whether our approach 


is useful or not is the degree to which the flare ribbons match the QSLs derived from unstable mod¬ 
els for erupting sigmoidal active regions. The validity of this approach is discussed in Section[3j 
The paper is organized as follows: In §2, we present the observations of the erupting sigmoidal 
active regions. In §3, we discuss the specifics of the magnetic modeling, and §4 details the QSL 
calculation. In §5, we present our main results about the association of QSLs with flare ribbons. 
Concluding remarks, along with a preview of Paper II, are given in §6. 


2. Observation 


In Savcheva et al.1 (2014), we presented a sigmoid catalog spanning the duration of the Hinode 
and SDO mission^]. For this paper, we have chosen seven sigmoidal active regions selected from 
the catalog. We use the recently-proposed naming convention for solar targets endorsed by the 


1 Hinode & SDO Sigmoid Catalog: http://aia.cfa.harvard.edu/sigmoid/ 
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Astrophysical Journal, Solar Physics , and the International Astronomical Union, under which our 
7 regions are labeled: 

SOL2007-02-12T07:30:00L075C297 

SOL2007-12-07T04:20:00L085C296 

SOL2010-04-08T02:35:00L110C176 

SOL2010-08-07T17:55:00L100C002 

SOL2012-05-08T09:21:00L108C232 

SOL2013-03-15T05:46:00L116C077 

SOL2013-04-llT06:55:00L110C071 



Regions were chosen based on the following criteria: 


1. Being eruptive sigmoidal regions, which are taken as a dir ect evidence for the existence of 
magnetic flux ropes before the eruptions ( Green et, al.1 2011). 


2. Having very good observational coverage around the flare time. 

3. Having well-defined double ribbons on both sides of the PIL, i.e. having two-ribbon flares. 

4. Having observations of other flare- and CME-related features, such as transient coronal holes 
(dimmings) and flare arcades. 


This sample does not include all of the catalog regions that fit our criteria. The case-by-case 
numerical modeling and topology analysis that we performed on each region (see Sections [3] &: [4]) 
represents a substantial and time consuming effort. While this sample of seven regions remains 
statistically limited, it far exceeds similar published analyses of flaring regions (cf. Sect. [I]). 

It is important to emphasize that we select regions based purely on the quality of the data and 
on the observational properties of the EUV and X-ray emission before and during the flare. None of 
our selection criteria are directly based on the magnetic field properties of the region. The selection 
does not presume how well the magnetic models fit the data or how well the ribbons match the 
QSLs. The topological analysis was performed blindly, using only the magnetic held model and 
without taking in to account the EUV and X-ray emission of the individual region. We present all 
regions in the original sample, some of which are not perfectly matched. The list of the regions 
by date, including the time ranges over which different flare and CME features were observed, are 
presented in TableHJ 

As can be seen from Table[lJ our sample contains both weak (B- and C-class) and some stronger 
flares (up to M6), so we can demonstrate the correspondence of QSLs and ribbons in different types 
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of magnetic environments. Both STEREO flares are B-class flares that do not display clear motion 
of the ribbons, probably due to the low intensity of the flare. Our sample does not contain any 
X-class events because all suitable X-class flares have excessive saturation that precludes a clear 
examination of the flare ribbon evolution. 


In this study, we are concerned with the initial moment when the full lengths of the flare 
ribbons appear. This moment is defined by eye by inspecting image sequences in several AIA or 
STEREO channels. The ribbons appear at different times and to different extents in the various 
channels since they are sensitive to different temperatures. For the SDO regions, we rely primarily 
on the 304 A channel since it is often less saturated and the ribbon shape is not confused by 
overlying loops as in the other channels. We did not make the traditional choice to use 1600 A 
images because not all flare ribbons were seen in their entirety in this channel and because the 
ribbons generally appeared later in the 1600 A channel (i.e. the 304 A channel captures hotter 
plasma than the UV channels), and we wanted to obtain a moment as close to the beginning of 
the flare as possible. Because of the more limited observation cadence (5 min) for the STEREO 
regions, we chose the channel that showed these particular ribbons the best, i.e. the 195 A channel. 
In Fig. [0 and El we display the flare ribbons in all regions at the times when the ribbons are seen 
first at their full extent. This time, tribbons (also listed in Tabled]), is the reference time used in 
Section[5] to make the overlays between the QSLs and flare ribbons. 


Some ribbons brighten more or less simultaneously over roughly their whole length and others 
grow in length as time progresses. The moment, tribbons-, that we use to compare the ribbons with 
the magnetic model is the first moment when the ribbons are ob served over their full e xtent. In 
class i cal two-ribbon flares, ribbo ns display two kinds of motions (lAsai et al.112002 . 120041 ; iFletcher 


2009; 


Qiul 12009 ; iQiu et al . 1201 01). One is a “zipper” motion (fast elongation) parallel to the PIL 
during the impulsive phase of the flare with velocities of 10-100 km s^ 1 . The other is an expanding 
(separating) motion perpendicular to the PIL during the gradual phase of the flare with velocities of 
1-20 km s -1 . At t r ibbons, most of the zipper motion has taken place, but the expanding (separating) 
motion has not started or is at a very preliminary stage. The ribbons in both of the STEREO events 
did not show significant motion over the duration of the flare, probably because of the weaker flare 
intensity. For these events, because of the lower cadence, it was also harder to strictly pinpoint 
tribbons- As can be noted from Tabled! t r ibbons is usually within 10 min after tcoES-, of the flare 
start time measured by GOES. 


In Fig. O we show pre-flare observations from XRT for each of the seven regions, revealing their 
sigmoidal nature. The XRT observations are displayed close to the time, (listed in Table[2]), 
at which the 3D magnetic field models have been produced for each region. Most images show 
S-shaped loops that outline the sigmoids. All sigmoids on the figure are long-lasting with the 
exception of SOL2013-04-11, which was transient and appeared a couple of hours before the GOES 
flare. For this region, the XRT image is not taken at a time useful for constraining the model, i.e. 
no S- or J-shaped loops are visible, so we use AIA 335 A instead. 
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As we will discuss in Section^ the models are constrained by a number of coronal loops traced 
from these images. The images were chosen to show coronal structures that did not show any 
flare-related changes yet but were also representative of the immediate pre-flare configuration. For 
example, some of the loops seen in the figure lifted up as the flare was starting and others faded 
away, giving place to the flare arcade. 


3. Magnetic Field Modeling 


Since these regions are sigmoidal, we can gain insight into their 3D magnetic field structure by 
building NLFFF models, which are the only models that can represent the core field of a magnetic 
flux rope and the overlying potential arcade (|Savcheva &; van Ballegooiienl 120091 ) . We assume that 
the coronal field is force free, hence that V x B = aB, where B is the magnetic field and a is the 
torsional parameter that is constant along a given field line. A force-free magnetic field is “non¬ 
linear” when the torsion parameter is a function of position. Different field lines have different 
values of a, and a is as close to a constant along the field line as possible. In our models, the 
cores of the sigmoids are represented by sheared and twisted fields, while the outlying fields have a 
nearly potential configuration. The relaxation method ensures that the resulting field is force-free 
and non-linear. 


We have used the flux rope insertion method (Iv an Balleg ooii en 20041 ) to produce models of 
the magnetic field in each of our observed cases for the instants, t mag , preceding the eruptions 
by about 0.5-1 hour (see TablcQ] and [2]). Briefly, the flux rope insertion method consists of the 
following steps: 1) A global potential field extrapolation is performed based on a low-resolution 
synoptic SOLIS magnetogram. 2) A region of interest is selected from a high-resolution MDI or 
HMI magnetogram. 3) A modified potential field extrapolation is performed in the high resolution 
region with side boundary conditions set by the global potential field extrapolation. 4) A flux 
rope is inserted following the path of an EUV filament seen STEREO 195A or AIA 304A images. 
5) A grid of models is created with different combinations of axial and poloidal fluxes. 6) The 
configurations are relaxed to a force-free state using magnetofrictional relaxation with some amount 
of hyper diffusion, which acts to smooth gradients in a. 7) Field lines from each model are matched 
to observed coronal loops selected from X-ray or EUV (94 A or 335 A ) images. We choose loops from 
the core of the active region where the field is strongly non-potential, usually representing the spine 
and elbows of the sigmoid, as well as any linear, sheared, and more potential loops, thus spanning 
the whole magnetic configuration as well as possible. Consequently, this method produces a 3D 
magnetic field struc ture that i s heav ily c onstrained by observati ons. The details of the method are 
further discussed in ISu et all ( 2011 1 and Savcheva et al. ( 2012a ). We have a quantitative measure 
of the goodness-of-fit for eac h model. The goodness of fit criteria are explained extensively in 
Savcheva Sz van Ballegooiien (120091 ). By looking at the range of solutions near the best-fit model, 
we can estimate the robustness of the best fit, or how sensitive the models are to the choice of 
initial parameters. 



















The instances for which we build the models correspond to the times just before the first flare- 
related morphological changes start to appear in the AIA/STEREO images. For all regions, we 
calcul ate a grid of magnetic models coveri ng a range of flux rope axia l and poloidal fluxes, as we 
did in ISavcheva fc van Ballegooiienl (J2009j) and ISavcheva et al.1 (1201 2d ) . We have chosen the grid 
to cover models with little flux (e.g. <&£«; = 1 x 10 20 Mx and F po i = 5 x 10 9 Mxcm” 1 ) that produce 
sheared arcades to very highly twisted (large amount of polidal flux, e.g. Fp 0 \ = 1 x 10 11 Mxcrn 1 ) 
and/or sheared (large amount of axial flux, e.g. <L ax j = 5 x 10 21 Mx in some cases) fields that 


exceed the stability limit for the particular region. Su et ah] ((2011) showed that the stability of the 
configuration is tightly controlled by the amount of axial flux and that the fit of the models to the 
images is weakly dependent on the poloidal flux. We issued best-fit models that closely match the 
observed loop configurations in XRT synoptic images and AIA 335A and 94A channels when XRT 
was not available. This shows that our models are robust enough to reproduce excellent fits to high 
resolution AIA observations in the cor e of activ e regio ns. The first such NLFFF model constrained 
by AIA observations was provided by ISu et al . (2011,1). 


The flux rope insertion method produces models that are on a spherical wedge grid. However, 
the QSL computation is performed in Cartesian coordinates, so we transform the spherical to 
Cartesian coordinates of the grid as we did in S12. 

Although in previous studies we were concerned with the best f i t models that represented the 


equil i brium coronal magnet ic field preceding flares (ISu et al 


2009 


2011 


Savelieva & van Ballegooijen 


Savelieva et ah 2012q), here (and in a forthcoming paper) we are more interested in under¬ 


standing how the field evolves from a static to a dynamic flaring state and in what the magnetic 
field structure is at the flare onset. The flux rope insertion method produces a stable relaxed model 
when the forces on the flux rope bal ance. If no b a lance is found, the solution continues to evolve as 
the number of iterations increases. Kliem et al.l ( 2013 1) demonstrated that the stability boundary 
identified in our NLFFF code is reproduced by a full 3D MHD solution. This gives us confidence 
that the marginally stable solutions found by our NLFFF method are reasonable representations 
of the magnetic field at the onset of the flare instability. 

As discussed earlier, if too much flux is used in the flu x rop e, a relaxed equilibrium state may 
never be reached by the magnetofrictional relaxation (jSu et al, [2011|) and the model is unstable. 
The flux rope in these u ns table models continues to expand and rise with continued iterations, as 
demonstrated in[Su et ajj ((2011). The unstable models are hence not force-free since since residual 
forces exist in the volume that drive the expansion of the flux rope. In most cases, the best-fit 
model based on the pre-flare configuration is unstable or marginally stable. In the marginally stable 
cases, we add a unit of axial flux (1-2 x 10 2 ° Mx) to the axial flux of the best-fit model, which pushes 
the models over the edge of stability. 


There are usually no suitable coronal loops that can be used to match the unstable model to 
the data. The flare l oops are not a go od choice for fitting because they undergo the strong-to-weak 
shear transition (jA ulanier et alJ 2012l), and the match will depend on the iteration number, which 
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is a central topic of Paper II (see for details). We choose the closest unstable model to the best-fit 
model and further relax it with no form of diffusion except a numerical one. If the model is stable 
with the current relaxation scheme we are using, it converges to a NLFFF state at about a relaxation 
iteration number of 30,000. We save every 10 000-15 000 iterations in order to track the evolution 
of the magnetic configuration with further relaxation. In the case of the unstable models, th e flux 
rope keeps expanding and forms an HFT underneath it (ITitovI 120071 : ISavcheva et al.ll2012al lbl). In 


Tabled we have given the times, input and output parameters of the unstable models we use for all 
regions. The input parameters are the axial flux (d> ax j) and the poloidal flux per unit length ( F po i ), 
which when multiplied by the length of the flux rope path results in the total poloidal flux given 
in the table, F po i^ ot . The output parameters are potential field energy ( E pot ), computed before the 
introduction of the flux rope, the free energy ( Ef ree ), computed by subtracting the total energy of 
the NL FFF field a nd E, nnf . The last column is the relative helicity, computed from equation (B2) 
from iBobra et al.l (120081 ). Sample field lines from the best fitting models are shown in Fig.U One 
can compare Fig. [3] and Fig. [4] to see that there is a very good qualitative match of the field lines 
to the X-ray emission. 

In principle, modeling an unstable eruptive configuration with NLFFF is not easy, even when 
using vector fields to extrapolate the coronal field. This is in part because the photospheric magnetic 
field changes only very subtly in respons e to the fla re, an d this change is mostly observed in the 


transverse component of the vector field ((Sun et_ah 


2012 


Wang et al.ll2012l L In our case, we use 


LoS magnetograms in which the change is virtually undetectable. During the relaxation process, 
the vertical field distribution B z at the photosphere stays fixed and equal to the pre-eruptive LoS 
magnetogram that we used to obtain the best fit model. The B x and B y at the photosphere stay 
equal to zero in the boundary layer during the relaxation, enforcing strict line tied conditions 
and preserving the initial purely vertical boundary condition. Even if the change in the LoS field 
during a flare was apparent, it would be most probable that we cannot choose the initial model and 
relaxation pa rameters precisely enough to represent this chang e in the final 3D field. There have 
been studies (jJiang et al .1120121 : iLiu et al .II2Q14J : IZhao et al.ll2014f ) of the quasi-static evolution of the 
magnetic field covering the time of a flare by utilizing vector magnetogra m extrapola t ions, which 
in principle can capture the flare-related changes in the configuration. In Zhao et ah |2014|) they 
show the existence of a flux rope before and its absence after the flare. There is little discussion 
of how the flare-related changes to the photospheric vector field affect the coronal field resulting 
from the extrapolation. This is a fundamental difference between the flux rope insertion method 
and the vector NLFFF extrapolations in that after having identified a stable or marginally stable 
solution that matches the observations, we can perturb that solution so it becomes unstable. 

Although the unstable models are not fully self-consistent, using a marginally unstable model 
allows for a description of the flaring magnet i c con figurations close to the expected real 3D magnetic 


field in each reg ions at t r jhi, nns . Kliem e_t al. ((2013) performed a data-constrained MHD simulation 


initiated by the Su et al.1 (|201ll ) NLFFF models as initial conditions and showed that the dynam¬ 
ics of unstable models in the magnetofrictional relaxation is close to the dynamics of an unstable 
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flux rope in the MHD simulation. The magnetofrictional relaxation, being far more economic in 
numerical resources than an MHD simulation, allowed us to study marginally unstable models for 
all seven studied regions, which would have been computationally prohibitive using an MHD simu¬ 
lation. The marginally unstable models therefore allow us to obtain an economic and satisfactory 
representation of the flaring active region field. 

It is important to emphasize that although possible, we do not use any additional free param¬ 
eters or constraints on the models based on the eruptive or post-eruptive coronal or magnetogram 
structures. All parameters and constraints are applied using the pre-flare configuration. We are 
using the predictive capabilities of the standard model to test the quality of our reconstruction 
method. The topological analysis of our models, which was performed independently of the model 
production, allows us to predict the location of the ribbons. If the ribbons do not match the QSLs 
at all, we cannot say anything about the eruptive structure of the region and the reconnection 
site. On the other hand, if there is a good qualitative match, we can be confident that we have 
found a solution that represents the observations well and supports the validity of our theoretical 
prediction. 


4. The QSL calculation 


As mentioned in the introd uction, based on previous work fe.g. lSavcheva et al.ll2012bl : iLiu et al 


2014 ; Janvier et al.ll2013 , 2014 ) and the standard flare model, we expect that the flare ribbons match 
the QSLs in the low corona. In order to determine the fits between the magnetic model and the 
observed ribbons, we per form a top ol ogical analysis of the 3D magnetic field data by computing 
the squashing factor (Q) (Ti tov et al . 2002 i in diff erent 2D cuts of the domain. To compute these 
maps, we use the method of P aria t & Demo ulin (120 1 21 ) that was successfully applied to a single 
NLFFF model in S12. This is an iterative method that can take any cut through the domain and 
first compute a map of the QSLs with the original resolution of the model, which is 1.5 x 10~ 3 -Rq, 
followed by increasing the resolution where Q is high and repeating the process until convergence 
to a similar value of Q is reached. 


For each unstable model of the studied regions we, calculate the Q-rnaps at four heights in the 
corona: 2100km, 3150 km, 4200 km, and 5250km (z = 2 — 5 in simulation cells). The reference 
boundary is at z = 2, since below this height the intrinsic complexity of the magnetic field is 
too large, i.e. a large number of bald patches that exist make the converging QSL computation 
very slow, as discussed in S12. Taking the reference boundary for the QSL computation as the 
photospheric layer only focuses on topological features that are situated very low and prevents 
clear localization of the topological structures that extend into the corona and are involved with 
the large-scale flare dynamics. Taking the reference field a few pixels above would correspond to 
a chromospheric-like region where the plasma beta would be closer to unity but where line tying 
would still be effective. This also allows us to pick out significant topological features in the quiet 
sun where the held is usually very week, highly fragmented, and containg interspersed polarities. 
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Once the Q factor has been computed from an unstable model based on magnetic held data 
obtained at t mag , we overlay each of these QSL maps with an image of the flare ribbons at t r ibbon- 
The overlay is performed in such a way that the Cartesian QSL map is projected onto a sphere 
and the height at which the QSL map is calculated is taken into account in the projection. We 
overlay each map for the four different heights onto the image by shifting the QSL map coordinates 
forward in time to match the time difference between the moment when the model was computed 
and the time of the flare ribbon image, which in some cases was up to two hours. In addition, we use 
prominent magnetic features, such as sunspots and large, peculiarly-shaped magnetic concentrations 
to further adjust the alignment. The co-alignment precision can vary between different location 
in each co-aligned image because of the transformation from the spherical data to the Cartesian 
one. Hence, for each studied region, the coalignment may be better in a portion of the image than 
another. Overall, we estimate that this process gives us an alignment precision of about 5-10 AIA 
pixels or 3 // -6 ,/ on average. 


Given the complexity of both the QSL maps and the flare ribbon observations, the existence 
of a large number of QSLs related to small scale features, and the presence of multiple secondary 
brightenings, we only focus on a qualitative comparison of the QSL maps with the ribbon mor¬ 
phology. We determine the quality of the fit by eye employing a few basic criteria. The QSLs 
should first capture the portion of the ribbons parallel to the PIL (the straight part of the J-shaped 
ribbon). This section is indeed directly related to the extent and height of the HFT located below 
the erupting flux rope. Its dynamics are related to the rise of the flux rope and the reconnection 
developing along in the CME current sheet. Good alignment between the QSLs and the flare 
ribbons indicates that the central axis of the twisted flux rope is located at the right height and 
follows the right shape. Flare ribbons are frequently observed to not be fully facing each other, with 
positions slightly shifted in the direction of the PIL. This shift is related to the shear present in the 
pre-eruptive configuration. If the QSLs are able to capture such shift in position, it indicates that 
the modeled field contains the right amount of shear around the flux rope. The second criterion 
to evaluate the quality of the fit is the hook of the J-shaped flare ribbons, the sections going away 
from the PIL. The position of the hook is related to the position of the centr al axis of the twisted 
flux rope (jPariat Sz Demoulinll2012l ; ISavcheva et al.ll2012bl : lJanvier et alJl2013l ). Its shape is related 
to the amount of twist present in the rope: a more twisted flux rope will present a more pronounced 
hook, eventually closing in on itself for very large twist. Finally, we hope to capture smaller scale 
detail in the shapes of the ribbons. 


In most cases the maps at z = 2 Mm were still too complex with many small-scale features, 
although in the case of SOL2007-12-07 and SOL2010-04-08 these were the best matching maps. 
Therefore, we usually chose the maps at z = 4 Mm as the best match (the reference boundary for 
the computation of Q always being at z = 2 Mm). These best fitting maps are shown in Fig. [5] 
in a logarithmic scale. Note that we reach values of 10 10-12 in Q in these maps, which makes 
us confident that we capture the important QSLs that bind the flux rope and quasi-separate it 
from the surrounding field, as shown in S12. In this figure the highest values of Q are in reddish 
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shades, the background Q is blue, and areas where one end of the held line was located outside 
the computational domain (so Q could not be computed) are white. In some cases, the QSLs are 
somewhat S-shaped (SOL2010-08-07) and most are 2J-shaped (SOL2007-02-12, SOL2010-04-08, 
SOL2012-05-08, SOL2013-03-15). 

The chara cteristic 3D shape of the QSLs associated with a twisted flux rope is presented in in 
Fig. [6] (see also Titov et al. 2002 : Savcheva et al. 2012b : Pariat &; Demoulin 2012 1. Fig. [6] presents 
a synthetic description of the QSL model of the SOL2007-02-12 sigmoid. We show the electric 
current density distribution in a horizontal cut at level 4 Mm and at the location of the blue line 
- a current cross-section through the flux rope. We use the current as proxy for the locations for 
QSLs since in S12, we showed there is a very close correspondence although the current is more 
diffuse than the QSLs. The pinching under the flux rope is apparent in the figure where the HFT 
is. It is evident that since this is an unstable model, the HFT is already at some height above the 
photosphere. In green we show the location of three horizontal cuts through the flux rope, and on 
the right the schematic view of the QSL shapes are given depending on where the cut is taken. 
From this it is evident that if the horizontal cut is taken at the HFT, the QSL will look S-shaped, 
and if the cut is below it, it will look 2J-shaped, i.e with the straight parts of the Js parallel to 
each other and the PIL, and the hooks facing each other at the far ends. All the maps that do not 
show S-shaped QSLs show 2J-shaped QSLs corresponding to the cut under the HFT. 

The correspondence between current layers and QSLs is further demonstrated in Fig.0 where 
we show the horizontal distribution at the height of the best-fit QSL maps with the QSLs overlaid. 
The QSLs are selected to be above the same threshold as in Fig. [8] but also to be in areas where 
the magnetic field is more than 15-30 G for the different regions. This filtering allows us to remove 
a large part of the QSLs that are created from small scale polarities. Compared with Fig.[5l we 
highlight only the QSLs that are going though intense magnetic polarities, or where energy and 
currents are likely to be stronger. It is clear that indeed most QSLs overlay the edges of the current 
concentrations. This is to be expected since, as it can be seen from Fig. [6] top panel, the domain 
under the HFT is filled with intense current bound by the legs of the HFT. As shown in S12 and 
in paper II, most flux ropes show a hollow core current distribution, meaning that the current is 
concentrated on the edge of the flux rope. Then, in the flux rope, these ridges in the current are 
overlaid by the flux rope-binding QSL that separates it from the overlaying field. However, it is 
important to note that QSL theory does not expect an exact correspondence between the maxima in 
the current concentrations and QSLs since the exact relationship between the current and QSLs is 
complicated, not ana l ytically determined, and d ependent on the exact motions in the photosphere 
( Aulanier et al. 2006 : Wilmot-Smith et al. 2009 1. 


As theoretically predicted, our magnetic field models show the most intense currents accumu¬ 
late where the connectivity gradients are strongest. We showed previously in S12 that the flux rope 
has a hollow core current distribution with current concentrated at the edge, as can also be seen 
from Fig.[6l This current concentration coincides with the QSL t hat w raps around the flux rope, 


which is the case for all regions that we study. In Janvier et ah (|2013f ). this correspondence was 
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demonstrated in the 3D MHD simulation and in Janvier et alj |2015) for data-constrained QSLs 
and observed current ribbons. 


5. Flare Ribbons and CME Topologies 

Several studies have shown the correspondence between photospheric QSLs and ribbons in case 
by case examples (see Sect. [[]). Here, for the first time, we show good-quality matches between QSLs 
and flare ribbons in seven different regions containing flux ropes that exhibit classical two-ribbon 
flares. 

In Fig. [8] and 0 we show the matches between the QSLs (from Fig(5]) and the flare ribbons 
(from Fig (I]and [2]). In the left columns, the overlays between the images and all QSLs with Q values 
above 10' 5 ~ 6 are displayed (the Q threshold is marked above each panel). In the right columns, we 
overlay the flare images with QSLs with Q values only in areas where the current is above a certain 
threshold (the j threshold is also marked above each panel). 

In the following, we will not discuss whether a given QSL is associated with a given flare 
ribbon. In the real observed regions, while a QSL may be present, no current may accumulate at 
these QSLs and thus no flaring activity will be observed. We will thus mainly discuss whether an 
observed ribbon is associated with a QSL or not. The use of the threshold on the electric currents 
allows us to display mainly the QSLs along which currents have accumulated in sheets (FigjTj). As 
discussed in Section 4, these QSLs are most likely to host magnetic reconnection and hence show 
flare ribbons. By using this selection criterion, we highlight the QSLs that our models predict 
to be the preferential sites for the appearance of flare ribbons. A posteriori, we do observe that 
the overlays with the current criterion (right panels of Fig(8] and (9|) give a better correspondence 
between the ribbons and the QSLs than with the magnetic strength criterion as in Fig(71 In the 
left panels of Fig(8] and [9l numerous high-Q regions are observed to be associated with no EUV 
emission. In the real observed regions, while a QSL was actually present, the dynamics did not lead 
to any current accumulation at these QSLs and thus no flaring activity was observed. The better 
correspondence obtained with the right panels is a further demonstration of the overall predictive 
capabilities of the marginally unstable modeling method. 

The height at which the best-matching QSL maps are taken vary from 2 = 2 — 4 model cells 
or 2100 km to 4200 km. The times at which the overlays are made are those shown in Fig. [T] and 
(3 i.e. the first instance, t r ibi , ons , at which the whole ribbons appear. The iteration at which the 
QSL maps are taken are between 30,000 and 60,000 for the unstable models that are closest to the 
best-fit model. We will show in Paper II that as the relaxation proceeds, an unstable model shows 
an elevating HFT with legs of that move apart in time (relaxation iteration) and map the motion 
of the expanding flare ribbons. 

Based on the criteria for the quality of the match between the ribbons and the QSLs, we 
now discuss each active region emphasizing the parts of the ribbons that are well captured as well 
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as the parts where the QSLs fail to follow the ribbons, along with some possible reasons why. 
Particularly, we look for a match in the shape of the ribbons and QSLs in the straight (parallel to 
the PIL) part as well as the hooks. The separation of the QSLs/ribbons depends on the height of 
the of the reconnection site and the geometry of the legs of the HFT. These control the separation 
of the footprints of the HFT at chromospheric heights, where the ribbons are produced. The 
observations do not give us any information about the height of the reconnection site at the onset 
of the flare. Hence, we simply choose an iteration from the magnetofrictional evolution of the flux 
rope that has the HFT at the appropriate height to match the initial ribbon separation from the 
observations. As we show in paper II, the subsequent ribbon separation in time is matched by the 
rising HFT in the models. The evolution of the HFT with iteration number depends mostly (if not 
solely) on the amount of axial flux. The poloidal flux, according to the standard 3D model, rather, 
controls the tightness of the hooks but not the separation of the ribbons in the part parallel to 
the PIL. For a future study, we have planned to explore in detail the influence of different initial 
parameters on the shape of the QSLs. In this study, however, we restrict ourselves only to showing 
the closest marginally stable or unstable model to the best-fit pre-emption model, which contains 
the predictive power of this method. 


SC)L2007-02-12T07:30:00L075C297 (Fig.El first ro w): This is the same region we s tudie d 
in S12, the time evolu tion of which was explored in ISavcheva fc van Ballegooiienl ([20091); 


Savcheva et al. (j2012d ). For this weak B-class flare, the QSLs do not match well the flare 


ribbons, in particular in the east part of the sigmoid. This is probably because the flux rope 
path that we insert is too long, since it could not be constrained by the STEREO images very 
well. The endpoints of the dark EUV filament were not well discernible. The weak strength 
of the flare also corresponds to a limited intensity of the flare brightennings. We again 
emphasize that in all cases we only use the pre-eruptive observational information to constrain 
the models, but in this case the match could be improved if one uses the round transient 
CHs, with center location around (x,y) =(470,-90) and (380, -160), to constrain the anchor 
points of the flux rope better. The presence of weaker magnetic fields compared to other 
studied regions also reduced the constraint on the magnetic field during the magnetofrictional 
relaxation. Fields lines are more free to connect in various places of the domain than in a 
region with stronger field. This region thus represents a case which is below the lower limit, 
in terms of flare class and magnetic field strength, to what one is able to model with a good 
cor r esp ondence. 


• SOL2007-12-07T04:20:00L085C296 (Fig. El second row): This is another relatively weak 
flare (B2) in a relatively quiet region. These ribbons are very well matched by the computed 
QSLs, which almost exactly overlay both the parallel parts of the J-ribbons around x = —270 
and the hooks around (-280, 0) for the Northern ribbon and x = —300 for the Southern 
ribbon. The y coordinate goes from -120 to 25. In addition, a QSL binds the region marked 
“TCH” for transient coronal hole in the corresponding image in Fig. El which is where the 
flux rope northern endpoint is situated and the dimming is observed as reduced emission. In 
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paper II, we will explore whether the dimming regions also match the QSLs outlining the 
footprints of the flux ropes. This analysis has large implications to extending the standard 
flare model to 3D in realistic coronal magnetic fields. 


SOL2010-04-08T02:35:00L110C176 (Fig.EJ third row): This is a B4 class flare occurring 
in a region with more intense magnetic fields. The central straight part of the ribbons, the 
bars of the two J’s around (-250, 440), are very well matched in this region. Although a nice 
eastern hook is visible in the QSL map around (-330, 375), the observed ribbons display a 
much narrower hook, i.e. the bar of the J and the straight part of the hook of the J are 
very close to each other around (-300, 400). This can clearly be seen in the movie provided 
in the online version. The hooks become more prominent at later times. Later, the eastern 
hook is relatively wide (but not as much as the one in the QSLs) and the western one is also 
visible. This region was extensively modeled by[Su et al. |2011), and the model is very well 
constrained. However, it is possible that the surrounding field limits the extension of the flare 
ribbons, while this behavior is not seen in the QSL. The dynamics of the QSLs in this region 
will be discussed extensively in Paper II. 


• SC)L2010-08-07T17:55:00L100C002 (Fig. [HI fourth row): This M class flare is also overall 
very well matched by the model. The ribbons in this region resemble more an S-shape rather 
than two Js. Both the north-east hook around (-580, 160), and the more parallel part around 
(-500, 130) of the ribbons are well represented with only slight deviations. The swirly ribbon 
coinciding with the sunspot at (-480, 70) in the south-west end of the main ribbon is not 
captured by the QSL map since field lines that originate in this area leave the computational 
domain and Q cannot be calculated for them; these areas appear white in the corresponding 
QSL map in Fig.[5j In the left panel that contains all QSLs, there are two QSLs that pass 
near the north-east hook, but after the threshold on the current is imposed (in the right 
panel) only one remains. These two QSLs spread later on and overlay the two spreading flare 
ribbons, as we will discuss in Paper II. 


• SC)L2012-05-08T09:21:00L108C232 (Fig.El first row): For this C class flare, the north¬ 
east hook around (-40, 375) as well as the straight part, centered at (10, 320), of the 2J ribbons 
are well matched by QSLs. However, there is basically no south-west ribbon hook where the 
QSL clearly forms a hook around (90, 240) at this early time. Later, if one refers to the movie 
in the online version one can see a hook forming at the location of the western hook in the 
QSLs. This region’s NLFFF model was constrained by an almost straight, slightly S-shaped 
filament observed in AIA 304A that clearly coincided with the S-shaped loops seen in XRT 
before the flare. However, as the eruption developed another branch of this filament lying to 
the east of the southern end of the first filament became apparent (with its Easternmost end 
at round (-60, 200)). A ribbon associated with the lift-off of this other filament appears in 
the images around x = —60 to x = 50 and y = 220, which is not associated with any QSLs. 
We performed alternative models, with different flux rope paths, but were unsuccessful at 
reproducing this part of the ribbon. The details in this modeling will be discussed in Paper 
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II. 

• SOL2013-03-15T05:46:00L116C077 (Fig.EJ second row): This M class flare displays a 
complicated pattern of ribbons with a lot of complexity in the basic shape of the ribbons 
and possibly more than two ribbons, but nonetheless a 2J-shape is clearly discernible. These 
2J-shaped ribbons are approximately well matched by the QSLs. The best match is found 
in the parallel part of the ribbons centered at (-180, 240). The hooks around (-230, 320) 
and (-60, 200) are more poorly matched due to their higher complexity, which is missed by 
the QSLs. The photospheric magnetic field distribution of this region, and SOL2013-04-11, is 
significantly more complex. These two regions also produced the strongest flares in our sample 
(M-class flares). These facts might be the reason for the complex shape of the ribbons. It is 
also possible that the overlying envelope of the sigmoid is not potential, but there is significant 
non-potentiality away from the flux rope, which can cause some of the mismatch. 

• SOL2013-04-llT06:55:00L110C071 (Fig.[9l third row): The situation with this region 
is very similar to SOL2013-03-15, i.e. strong flare, complex magnetic field distribution, and 
possible non-potentiality away from the main flux rope. However, in this case both the 
straight part (centered at (-240, 270)) and the hooks (around (-300, 290) and (-190, 220)) 
are very well matched, including the western hook which wraps around the sunspot at (190, 
235). The extent of the ribbons is well matched as well as the little curvy protrusion in the 
eastern ribbon at (-250, 290). 

Overall, the QSLs match the flare ribbons very well in six cases out of seven. This is a success 
for our approach since we have to deal with large complexity in the data-constrained QSL maps 
and the observed flare ribbons. In addition, these models are only constrained by the observed 
pre-flare loops and no information about the eruptive features such as ribbons, transient coronal 
holes, or flare loops, are used, although this would likely improve the match in some cases as we 
discussed above. 

In most cases, we find that the straight part of all ribbons is very well matched by the QSLs. 
Both hooks match well in four of the regions, and all regions have one hook matching. The hooks are 
generally harder to reproduce because this requires the most precise determination of the locations 
of the footprints of the filament that we use to guide our flux rope insertion. In addition, the shape 
of the hooks is dependent on the twist of the flux rope, as discussed above, and our models are 
less sensitive to the initial amount of poloidal flux in the flux rope, so the twist is less accurately 
determined. 

In the case where we do not have a very good match between the ribbons and QSLs, the model 
is not well constrained by the data before the eruption. The flare in this region is also the weakest 
and appears in highly dispersed magnetic polarities. 

Since we have used only pre-eruptive information about the coronal structure and surface 
magnetic held, the above, largely successful, comparisons between ribbons and QSLs show the 
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predictive power of our method. We can infer the eruptive topology and related features with a 
good success based only on the pre-eruptive parameters. This link will become even stronger in 
paper II, where we show that this way we can also match the observed evolution of the ribbons 
with time and the initial configuration and the strong-to-weak shear transition in the flare loops. 


6. Discussion and Conclusions 


In this paper, we showed the first comprehensive study of a selection of seven flaring sigmoidal 
active regions (SOL2007-02-12, SOL2007-12-07, SOL2010-04-08, SOL2010-08-07, SOL2012-05-08, 
SOL2013-03-15, and SOL2013-04-11) that display flare ribbons, transient coronal holes, and flare 
loops during solar flares. We select the regions purely based on the observations, i.e having 2 
J-shaped flare ribbons that appear in sigmoidal active regions and have a very good coverage of 
the pre-flare and flare features. Here we put the ribbons in a topological context, utilizing quasi- 
separatrix layer (QSL) maps obtained based on data-constrained magnetic field models. In order to 
compute the topology in each of the seven erupting sigmoidal regions, we use the flux rope insertion 
method to produce marginally stable and unstable magnetic field models to derive the 3D magnetic 
field structure at the flare onset. The use of unstable models in a magnitofrictional approach is 
justified by the fact that when used as initial conditions in an MHD simulation these unstable 
models evolve in a very simi lar manner a s in the subsequent iterations of the magnetofrictional 
relaxation process (see Kliem et alj 2013, for details). In this sense we use only observational 
information about the pre-eruptive magnetic field configuration and to predict the location, shape 
and extent o f flare- and CME-relate d fea tures. The QSL computa tion is performed in the same 


manner as m 


Savcheva et al.1 ( 201211 ) and lPariat &; Demoulinl (2012|). 


We show an very good match between QSL maps taken at low heights in the corona and 
chromosphere and the locations, extent, and shape of flare ribbons at the beginning of the flares 
in six out of seven cases that cover a wide range of flare intensities. Generally, the straight part 
of the J-shaped ribbons that are parallel to the PIL are better matched in all cases. To match 
the hooks, one requires very detailed knowledge of the location of the footpoints of the flux rope 
and the amount of twist contained in it. In this sense, our models are successful at reproducing 
both hooks in four of the cases. S012007-02-12 does not display a very good match in the eastern 
hook due to the poorly constrained path of the pre-eruptive flux rope. This study presents the 
largest sample of correspondence of active region flare ribbons and topological features derived from 
data-constrained magnetic field models. 


Being confident that our models match the observed flare features gives us the ability to 
analyze the flaring configurations. We can go further by using the “requirement” that the flare 
ribbons match the QSLs and apply it as a constraint on our models in order to understand the 3D 
reconnection process. For example, one can determine the height of the HFT at different locations 
along the flux r ope and the location of the flare. One can perform torus instability diagnostics, as 
we did in ( 1 Savcheva et al . 2012b, , S12), knowing the shape, size, and center of the flux rope from 
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the QSL cross-sections. One can look at the connectivities of field lines and their change during the 
eruptions. We reserve most of this analysis for Paper II, where we will perform extensive dynamics 
analysis on several consecutive iterations of the relaxation process and show the evolution of the 
QSLs and field line connectivities with iteration number (time). The analysis given in this paper 
and in Paper II has potentially very strong implications for the study of 3D reconnection in realistic 
coronal magnetic fields, which depart from the usual symmetry that many reconnection simulations 
and laboratory experiment employ. 


Previous studies have used LFFF and p oten t ial models to m atch observe d flare features to 
separatrices and QSLs (see Sect.[[]). Recently [Liu et ah (201411 and IZhao et al.1 (2014) used stable 
NLFFF models of complex regions with flux ropes. However, they did not interpret the connec¬ 
tion between the flare ribbons and low-lying QSLs in the same way as we do here. While such 
NLFFF methods are useful for studying magnetically complex regions that do not show a clear 
flux rope evidence, th ey are largely limited in their applicability (flux rope systems) as discussed in 


Schrijver et al. (2006), and they are intrinsically tied to the pre-flare configuration, since the surface 


magnetograms do not usually change significantly in response to the flare. On the other hand, by 
using unstable flux rope configuration in a magnetofrictional method we have the ability to follow 
the rise of the flux rope as the flare (reconnection) progresses, although we caution that this is not 
a self-consistent approach. We simply use it as a simpler model that maps the 3D magnetic field 
structure of the observed phenomena surprisingly well. To be fully consistent, one needs to use a 
data-driven MHD simulation or a time dependent magnetofrictional approach, which is planned for 
a future study to confirm this result. 


Most of our knowledge about reconnection comes from effectively 2D reconnection exper i¬ 
ments, and only recently efforts to understand 3D magnetic recon nection have been pursued (IPriest 


2011 


Pontin 


2011 


Shep herd fe C assak 


2012 ; 


Janvier et al.l 2014), revealing a much wider range of 


dynamics. Very little is known about the actual properties of the solar plasma during the recon¬ 
nection process on the Sun. Many authors have studied the global and local reconnection rates 
and timescales of s olar reconnection during s olar flares based on analysis of t he ribbon motions in 


the chromosphere (IFletcher fc Hudson 


Qiu et al, _2002; S ab a et al 


200 ' 


3 


2001 


Isobe et al 


Xie et al. 200HQiu et al.l l2010l : iMiklenic et al.l 120071 : 
200,4 Lling et al.lljOOrJ ) Generally, they use measure¬ 


ments of the magnetic flux at the photosphere, the velocity of the ribbons parallel to the PIL, 
and the 2D approximation of the standard flare model to determine the reconnected flux and the 
electric field and Poynting flux at the reconnection site. In a future paper, we plan to test the va¬ 
lidity of this approximation. However, these studies are purely based on the ribbon dynamics and 
are not able to capture the dynamics of the 3D coronal magnetic field, which is actually involved 
in the reconnection process. For example, to calculate the energy release rate from reconnection, 
one needs to combine observations with knowledge of the location and size (i.e. length, the width 
is generally under the resolution) of the reconnection current sheet, which can only be provided 
by data-constrained magnetic field models and extrapolations or by data-driven MHD sim ulati ons 


with the appropriate NLFFF initial conditions (and is still only an approximation). Aulanier et al 


































































-19- 


|2000, 2012) and ISchriiver et ah ([201 lj) utilize idealized MHD simulations to interpret the appear¬ 
ance of the flare ribbons. Although such studies show a qualitative similarity of the QSLs and flare 
ribbons and are valuable for determining the basic topology of the region, no actual estimates of the 
reconnection parameters can be obtained, which could be improved by the use of data-constrained 
and data-driven magnetic field modeling. 

So far, dynamics studies of the reconnection process have been either purely observational 
(reconnection rates, reconnected flux, etc.) or based on idealized MHD simulations. A quasi- 
dynam ic (as in the current approach), dynamic magnetofrictional (from sequence of magnetograms 


as m 


Kliemetal 


Gibb et al.l (201411. or a data-driven MHD simulation with a NLFFF initial condition (as in 


2013) will prove to be more valuable in the study of the evolving ribbons as the flare 


progresses and ultimately in the study of the process and properties of observed 3D reconnection 
on the Sun. 
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Table 1: Summary of the timing of observed flare features for all seven regions. 


ID 

Event 

GOES class 

tGOES 

Ribbons 

AIA 304A / 
STEREO 195A 

^ribbons 

Transient CHs 
AIA 193A / 
STEREO 195A 

Flare loops 
AIA 171 A/ 
STEREO 195A 

SOL2007-02-12 

12 Feb 2007 

<B 

07:40 

07:35-11:55 

07:35 

07:45-11:55 

09:05-10:45 

SOL2007-12-07 

07 Dec 2007 

B2 

04:20 

04:35-05:15 

04:36 

04:25-05:15 

04:55-06:45 

SOL2010-04-08 

08 Apr 2010 

B4 

02:35 

02:20-06:00 

02:39 

03:02-06:00 

02:49-06:00 

SOL2010-08-07 

07 Aug 2010 

Ml 

17:55 

17:52-20:27 

18:05 

18:12-19:02 

18:36-20:31 

SOL2012-05-08 

08 May 2012 

Cl.7 

09:21 

09:22-13:00 

09:30 

09:41-10:16 

09:56-12:58 

SOL2013-03-15 

15 Mar 2013 

Ml.l 

05:46 

05:57-09:00 

06:08 

06:11-gap 

gap-09:00 

SOL2013-04-11 

11 Apr 2013 

M6.5 

06:55 

06:27-08:00 

06:57 

07:00-08:00 

07:24-08:00 
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Table 2: Best-fit model parameters for all modeled regions. 


ID 

Event date 

tmag 

®axi 

[10 20 Mx] 

F pol 

[10 10 Mx cm -1 ] 

Epol ,tot 
[10 20 Mx] 

Epot 
[10 31 erg] 

Efree 
[1CT 1 erg] 

^rel 

[10 42 Mx 2 ] 

SOL2007-02-12 

12 Feb 2007 

06:41 

5 

5 

3.5 

3.64 

1.47 

1.08 

SOL2007-12-07 

07 Dec 2007 

03:21 

7 

5 

1.1 

1.28 

0.59 

0.34 

SOL2010-04-08 

08 Apr 2010 

02:00 

6 

1 

4.1 

5.29 

1.82 

1.85 

SOL2010-08-07 

07 Aug 2010 

17:00 

15 

-0.5 

-3.6 

29.7 

8.22 

-6.88 

SOL2012-05-08 

08 May 2012 

05:38 

7 

-1 

-4.0 

10.5 

2.83 

-2.53 

SOL2013-03-15 

15 Mar 2013 

05:00 

5 

-0.5 

-2.4 

53.8 

7.34 

-6.61 

SOL2013-04-11 

11 Apr 2013 

04:50 

7 

-0.5 

-2.5 

41.2 

9.26 

-4.78 
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Fig. 1.— STEREO (when AIA is not available) and AIA images for four of the regions showing the 
flare ribbons at t r «bbon> he. the first moment when they appear (see text about how we determine 
this moment). These images are used to match to the QLSs in Fig. [HI 
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Fig. 2.— AIA images for all regions showing the flare ribbons at t r iM>on- These images are used to 
match to the QLSs in Fig.[9j 
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Fig. 3.— X-ray images from XRT used to fit the NLFFF models for all regions. Notice the S-shaped 
and 2J-shaped loops in most regions that give the sigmoidal shape. 
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Fig. 4.— Some sample field lines traced from the best-fitting unstable models for all regions studied 
representing the 3D magnetic field structure. The red and green contours mark the positive and 
negative flux distributions in the photosphere respectively. The loops are overlaid onto the XRT 
images used to fit the models (see Fig. [3]). 
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Fig. 5.— QSL maps for all regions studied here. The QSL maps are in logarithmic scale in Q. 
The axes are given in model unites, in which one unit in x or y is 1 Mm. This corresponds roughly 
1.5" for the purpose of comparing with the observations figures. Notice the J and S-shaped QSLs 
associated with the flux ropes. 
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Fig. 6.— A horizontal map showing the current in the flux rope (lower panel) and a cross-section 
through the flux rope at the location of the blue line in the lower panel (upper panel). The pinching 
under the flux rope at the HFT is visible. The synthetic shape of the QSLs is shown in the right 
depending on where the horizontal cut is taken with respect to the HFT (shown schematically on 
the vertical cross-section through the flux rope). 
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Fig. 7.— Overlays between the vertical current density distribution in greyscale at the height of 
the best-fit QSL maps overlaid with QSLs with Q > 10 6 (red lines). The QSL are only plotted in 
areas with LoS magnetic field larger than 15-30 G. 
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Fig. 8.— Overlays between the flare ribbons (underlying image) shortly after the beginning of the 
flare and the QSL maps (red or blue curves) for the studied regions. All regions with Qthresh > 10 6 
are shown in the left column. In the right column we have plotted only the QSLs which are above 
the same values of Q but we have also required that the current density is above some value (usually 

jthresh 0.2). 
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Fig. 9.— Same as Fig. [8] for the rest of the regions. 






